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Abstract 

The mechanism which discriminates the pattern classes at the same A, is found. It is closely 
related to the structure of the rule table and expressed by the numbers of the rules which break the 
strings of the quiescent states. It is shown that for the N-neighbor and K-state cellular automata, 
the class I, class II, class III and class IV patterns coexist at least in the range, < A < 1 — j^. The 
mechanism is studied quantitatively by introducing a new parameter F, which we call quiescent 
string dominance parameter. It is taken to be orthogonal to A. Using the parameter F and A, 
the rule tables of one dimensional 5-neighbor and 4-state cellular automata are classified. The 
distribution of the four pattern classes in (A,F) plane shows that the rule tables of class III pattern 
class are distributed in larger F region, while those of class II and class I pattern classes are found 
in the smaller F region and the class IV behaviors are observed in the overlap region between 
them. These distributions are almost independent of A at least in the range 0.25 < A < 0.75, 
namely the overlapping region in F, where the class III and class II patterns coexist, has quite 
gentle A dependence in this A region. Therefore the relation between the pattern classes and the A 
parameter is not observed. 
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I. INTRODUCTION 



Cellular automata (CA) has been one of the most studied fields in the research of complex 

n 

systems. Various patterns has been generated by choosing the rule tables. Wolfram [1J has 
classified these patterns into four rough categories: class I (homogeneous), class II (peri- 
odic), class III (chaos) and class IV (edge of chaos). The class IV patterns have been the 
most interesting target for the study of CA, because it provides us with an example of the 
self-organization in a simple system and it is argued that the possibility of computation is 
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realized by the complexity at the edge of chaos 

Much more detailed classifications of CA, have been carried out mainly for the elemen- 
tary cellular automata (3-neighbor and 2-state CA) jj, [| , in which the patterns are studied 
quite accurately for each rule table. And the classification of the rule tables are studied 
by introducing some parameters □ □ □ . However, in this case, there are some confusions 
assification of the class IV CA. The result on the so to speak p — 1/2 task from 
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3| have been different from that of Mitchell, Crutchfield and Hraber[2|. In this case, 
the number of the independent rule tables are so small to treat them statistically and the 
symmetry of the interchange of the states "0" and "1" make the classifications of the the 
pattern classes more delicate than those of other CA with K > 3. Therefore, it may be 
worth to start with other models, in order to find the general properties of the CA. 

However the number of rule tables in N-neighbor and K-state cellular automata CA(N,K) 
grows like K K . Therefore except for a few smallest combinations of the N and K, the num- 
bers of the rule tables become so large that studies of the CA dynamics for all rule tables 
are impossible even with the fastest supercomputers. Therefore it is important to find a set 
of parameters by which the pattern classes could be classified , even if it is a qualitative one. 

Langton has introduced A parameter and argued that as A increases the pattern class 
changes from class I to class II and then to class III. And class IV behavior is observed 
between Cass II and Cass III pattern cto«flU Q. Tne A parameter represents rough 
behavior of CA in the rule table space, but finally does not sufficiently classify the quan- 
titative behavior of CA. It is well known that different pattern classes coexist at the same 
A. Which of these pattern classes is chosen, depends on the random numbers in generating 
the rule tables. The reason or mechanism for this is not yet known; we have no way to 
control the pattern classes at fixed A. And the transitions from a periodic to chaotic pattern 
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classes are observed in a rather wide range of A. In Ref. a schematic phase-diagram 
was sketched. However a vertical axis was not specified. Therefore, it has been thought 
that new parameters are necessary to arrive at more quantitative understandings of the rule 
table space of the CA Q- 

In this article, we will clarify the mechanism which discriminate the class I, class II, class 
III and class IV pattern classes at fixed A. It is closely related to the structure of the rule 
tables; numbers of rules which breaks strings of quiescent states. It is studied quantitatively 
by introducing a new parameter F, which we will call quiescent string dominance parame- 
ter. It is taken to be orthogonal to A. In the region 1/K < A < 1 — 1/K, the maximum 
of F corresponds to class III rule tables while minimum of F, to class II or class I rule 
tables. Therefore the transition of the pattern classes takes place somewhere between these 
two limits without fail. By the determination of the region of F, where the change of the 
pattern classes takes place, we could obtain the phase diagram in (A,F) plane, and classify 
the rule table space. 

The determination of the phase diagram is carried out for CA(5,4). It is found that the 
rule tables are not separated by a sharp boundaries but they are represented by probability 
densities. Therefore we define the equilibrium points of two phases where the two proba- 
bility densities of the pattern classes become equal, and define the transition region where 
the probability densities of the both pattern classes are not too much different from each 
other. By using the equilibrium points and the transition region, we draw a phase diagram 
in (A,F) plane. 

It is found that A dependences of the equilibrium points and transition region are very 
gentle, and they continued to be found at least over the range 0.125 < A < 0.75. It means 
that all the four pattern classes do coexist over the wide range in A. Our results for the 
distributions of these pattern classes in (A, F) plane do not support the well known relation 
between the pattern classes and the A parameter proposed in the Ref.J^j]. It will be shown 
that the results there, are due to the methods to generate the rule tables with probability 
A. 

In section II, we briefly summarize our notations and present a key discovery, which leads 
us to the understanding of the relation between the structure of the rule table and pattern 
classes. It strongly suggested that the rules which break strings of the quiescent states play 
an important role for the pattern classes. 
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In section III, the rule tables are classified according to the destruction and construction 
of strings of the quiescent states and we find the method to change the chaotic pattern 
class into periodic one and vice versa while keeping A fixed. We will show that by us- 
ing the method, the change of the pattern classes takes place without fail, in the region 
l/K < A < 1 - l/K. 

In section IV, the result obtained in section III is studied quantitatively by introducing 
a new quiescent string dominance parameter F. It is determined by using the distribution 
of class IV rule tables, which we call optimal F parameter. 

In section V, using F and A, we classify the rule tables of CA(5,4) in (A,F) plane. It will 
be shown that all the four pattern classes do coexist in wide range in A, contrary to the 
result of Ref. |3J] by Langton. The reason why he obtained his result will be discussed. 

In section VI, the rule tables of CA(5,4) are classified in the (A,F) plane, which provides 
us with the phase diagram. 

Section VII is devoted to discussions and conclusions, where a possibility of the transmis- 
sion of the initial state information and the classification of the rule tables by the another 
intuitive F parameter, will be discussed. 



II. SUMMARY OF CELLULAR AUTOMATA AND A KEY DISCOVERY 
A. Summary of cellular automata 

In order to make our arguments concrete, we focus mainly on the one-dimensional 5- 
neighbor and 4-state CA (CA(5,4)) in the following, because this is a model in which Langton 
had argued the classification of CA by the A parameter. However the qualitative conclu- 
sions in this article, hold true for general CA(N,K). These points will be discussed in the 
subsections III B. 

We will briefly summarize our notation of CA 

HQ. 

In our study, the site consists of 
150 cells having the periodic boundary condition. The states are denoted as s(t,i). The t 
represents the time step which takes an integer value, and the i represents the position of 
cells which range from to 149. The s(t, i) takes values 0, 1, 2, and 3, and the state is taken 
to be the quiescent state. The set of the states s(t,i) at a time t is called the configuration. 
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The configuration at time t + 1 is determined by that of time t by using following local 
relation, 



s(t + 1, i) = T(s(t, i-2), s(t, i - 1), s(t, i), s{t, i + 1), s(t, % + 2)). (1) 
The set of the mappings 

T(/i, u, K, p, a) = r], (/i, v, etc. = 0, 1, 2, 3) (2) 

is called the rule table. The rule table consists of 4 5 mappings, which is selected from a total 
of 4 1024 mappings. 

The A parameter is defined asj^. 

where Nh is the number in which r\ in Eq. (J2J) is not equal to 0. In other words the A is 
the probability that the rules do not select the quiescent state in next time step. Until 
section III, we set the rule tables randomly with the probability A. Our method is to choose 
1024 — Nh rules randomly, and set rj = in the right hand side of Eq. (J2J), and for the rest 
of the Nh rules, the rj picks up the number 1,2 and 3 randomly. The initial configurations 
are also set randomly. 

The time sequence of the configurations is called a pattern. The patterns are classified 
roughly into four classes established by Wolfram|lJ. It is widely accepted that pattern classes 
are classified by the A^; as the A increases the most frequently generated patterns change 
from homogeneous (class I) to periodic (class II) and then to chaotic (class III), and at the 
region between class II and class III, the edge of chaos (class IV) is located. 



B. Correlation between pattern classes and rules which break strings of quiescent 
states 

In order to find the reason why the different pattern classes are generated at the same 
A, we have started to collect rule tables of different pattern classes, and tried to find the 
differences between them. We fix at Nh = 450 (A = 0.44), because we find empirically 
that around this point the chaotic, edge of chaos, and periodic patterns are generated with 
similar ratio. By changing the random number, we have gathered a few tens of the rule 
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tables and classified them into chaotic, edge of chaos, and periodic ones. 

In this article, a pattern is considered the edge of chaos (class IV) when its transient 
length is longer than 3000 time steps. 

First, we study whether or not the pattern classes are sensitive to the initial configura- 
tions. We fix the rule tables and change the initial configurations randomly. For most of the 
rule tables, the details of the patterns depend on the initial configurations, but the pattern 
classes are not changed|l|. The exceptions will be discussed in the subsection VII A, in 
connection with the transmission of the initial state information. Thus the differences of the 
pattern classes are due to the differences in the rule tables, and the target of our inquiry 
has to do with the differences between them. 

For a little while, we do not impose a quiescent condition (QC), T(0, 0,0, 0,0) = 0, be- 
cause without this condition, the structure of the rule table becomes more transparent. This 
point will further be discussed in section IV. 

After some trial and error, we have found a strong correlation between the pattern classes 
and the QC. The probability of the rule table, which satisfies the QC is much larger in 
class II patterns than that in the class III patterns. This correlation suggest that the rule 
T(0, 0, 0, 0, 0) = h, h y£ 0, which breaks the string of the quiescent states with length 5, 
pushes the pattern toward chaos. We anticipate that the similar situation will hold for the 
length 4 strings of quiescent states. 

We go back to the usual definitions of CA; in the following we discuss CA under QC, 
T(0, 0, 0, 0, 0) = 0. We study the correlation between the number of the rules which breaks 
the length four strings of the quiescent states. These rules are given by, 

7X0,0,0,0,*) = /*, 

(4) 

T(z,0,0,0,0) = h, {i,h= 1,2,3). 
They will also push the pattern toward chaos. Similar properties of the rule tables had 



been noticed 
expandability 
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Wolfram and Suzudo with the arguments of the unbounded growth [l| and 



We denote the total number of rules of Eq. (jlj) in a rule table as A 4 . In order to study the 
correlation between the pattern classes and the number N±, we have collected 30 rule tables 
and grouped them by the number N4. We have 4 rule tables with ./V4 > 4 , 13 rule tables 
with N A = 3, 9 rule tables with N A = 2 and 4 rule tables with iV 4 < 1. When N4 > 4, all 
rule tables generate chaotic patterns, while when < 1, only periodic ones are generated. 




(e) (f) (a) (h) 



FIG. 1: The pattern classes at = 450. The quiescent state is shown by white dot, while other 
states are indicated by black points. Fig^a) corresponds to iV 4 = 4, Fig. ^b), Fig. [He), an d 
Fig. DJd), to N 4 = 3, Fig. [He), Fig. Q]T), and Fig. DJg), to iV 4 = 2, and Fig. DJh) corresponds to 
Ni = l, respectively. 

At iV 4 = 3 and N4 — 2, chaotic, edge of chaos, and periodic patterns coexist. Examples are 
shown in the Fig. ^ The coexistence of three pattern classes at iV 4 = 3 is seen in Fig. ^b), 
Fig. EKc) and Fig. DJd) and that of iV 4 = 2 is exhibited in Fig. QJe), Fig. [iff) and Fig. Dig). 

As anticipated, the strong correlation between A 4 and the pattern classes has been 
observed in this case too. These discoveries have provided us with a key hint leading us to 
the hypothesis that the rules, which break strings of the quiescent states, will play a major 
role for the pattern classes. 
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III. STRUCTURE OF RULE TABLE AND PATTERN CLASSES 
A. Structure of rule table and replacement experiment 

In order to test the hypothesis of the previous section, we classify the rules into four 
groups according to the operation on strings of the quiescent states. In the following, Greek 
characters in the rules represent groups 0,1,2,3 while Roman, represent groups 1,2,3. 
Group 1: T(/i, u, 0, p, a) = h. 

The rules in this group break strings of the quiescent states. 

Group 2: T(/i, u, 0, p, a) = 0. 

The rules of this group conserve them. 

Group 3: T(/i, u, i, p, a) = 0. 

The rules of this group develop them. 

Group 4: T(/i, u, i, p, a) = I. 

The rules in this group do not affect string of quiescent states in next time step. 

Let us denote the number of the group 1 rules in a rule table as N(gl). Similarly for the 
number of other groups. They satisfy the following sum rules, when is fixed. 

iV(<7l)+iV(<72) = 256, 

iV(<73)+iV(<74) = 768, . . 

N(g2) + N(g3) = 1024 - N h , 

N(gl)+N(g4)=N h . 

In the methods of generating the rule tables randomly using Nh (A), these numbers are 
determined mainly by the probability A; namely N(gl) ~ 256A, N(g2) ~ 256(1 — A) N(g3) — 
768(1 — A) and N(g4) ~ 768A respectively. Therefore they suffer from fluctuation due to 
randomness. 

The group 1 rules are further classified into five types according to the length of string 
of quiescent states, which they break. These are shown in Table H] The rule D5 is always 
excluded from rule tables by the quiescent condition. 

Our hypothesis presented at the end of the section II is expressed more quantitatively as 
follows; the numbers of the D4, D3, D2, and Dl rules shown in Table [Qwill mainly determine 
the pattern classes. 

In order to test this hypothesis, we artificially change the numbers of the rules in Table H] 



8 



TABLE I: The classification of the rules in group 1 into five types, where h ^ 0. 
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T(M,i, 0,0,0) =h 


12 






T(fx,j,0,0,m)=h 


36 


D2 


RP2,RC2 


T(i, 0,0, 1, a) = h 


36 






T(fi,j,0,l,(T) = h 


144 


Dl 


RP1,RC1 



while keeping the Nh (A) fixed. For D4 rules, we carry out the replacements defined by the 
following equations, 

T(0, 0, 0,0,i) = h-> T(0, 0, 0, 0, i) = 0, 
or T(i, 0, 0, 0, 0) = h -> 0, 0, 0, 0) = 0, 
^Xp, i/, j, p, a) = -> T(p, i/, j, p, a) = /, 



where except for h, the groups p, i/, p, a, j and / are selected randomly. Similarly the 
replacements are generalized for D3, D2, and Dl rules , which are denoted as RP4 to RP1 
in Table |l] They change the rules of group 1 to that of group 2 together with group 3 to 
group 4 and are expected to push the rule table toward the periodic direction. 
The reverse replacements for D4 are 

T(0, 0, 0, 0, i) = -> T(0, 0, 0, 0, i) = h, 
or T(i, 0, 0, 0, 0) = -> T(i, 0, 0, 0, 0) = h, 
T(fi, v, j, p, <t) = I -> T(p, v, j, p, a) = 0, 



which will push the rule table toward the chaotic direction. In this case, the groups h, p, u, 
j, p, and a are selected randomly. Similarly we introduce the replacements for D3, D2, and 
Dl, which will be called RC4 to RC1 in the following. 

By the replacement of RP4 to RP1 or RC4 to RC1, we change the numbers of the rules 
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FIG. 2: An example of the replacement experiments at = 615. Fig. Ufa) is obtained randomly 
with Nh = 615 (A = 0.6), by the method explained in subsection II A. Fig. Hfb) to Fig. Efh) are 
obtained by the replacements of the rule table of Fig. I2fa), which are summarized in Table UTI 

in Table U while keeping the A^ fixed. We denote these numbers A" 4 , A3, N 2 , and N% for 
D4, D3, D2, and Dl rules, respectively. By applying these replacements, we could study the 
rule tables which are difficult to obtain using only Nh (A). 

Examples of the replacement experiments are shown in Fig. |21 In this replacements, 
the RP4s are always carried out first, after that RP3s are done. The rule table of Fig. Ufa) 
is obtained randomly with A^ = 615 (A = 0.6). The numbers of the D4, D3, D2 and Dl 
are shown in line Fig. |2Ja) of Table |HI At Nh = 615 most of the randomly obtained rule 
tables generate chaotic patterns. We start to make RP4 tree times, then number of the 
D4 becomes A4 = 0. At this stage, the rule table still generate chaotic patterns. Then 
we proceed to carry out RP3. The chaotic patterns continues from N 3 = 22 to A3 = 15, 
and when N3 becomes 14, the pattern changes to edge of chaos behavior, which is shown 
in the Fig. Efb). Fig. |2fc) is obtained by one more RP3 replacements for Fig. Efb) rule 
table. It shows a periodic pattern with a rather long transient length. The pattern with 
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TABLE II: The numbers of the group 1 rules and numbers of the replacements for the rule table 
of Fig. Ufa) to change the patten classes shown in Fig. H 
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one more replacement of RP3 for the Fig. H c )> * s shown in Fig. Hd), where the transient 
length becomes shorter. These numbers of D4 rule (A/4) and D3 rule (A/3), and numbers of 
the replacements RP4s and RP3s for the Fig. H a ) rule table, are summarized in Table ILT1 

Similarly replacement experiments in which the RP4s are stopped at A/4 = 1 and A/4 = 2 
are shown in Fig. H e )> Fig. HO 5 an d Fig. Hg)? Fig. H n ), respectively and the numbers of 
group 1 rules and the replacements are also summarized in the corresponding lines in the 
Table ITT1 We should like to notice that the transitions to class III to class IV pattern classes 
take place at (N 4 = 0, N 3 = 13), (N 4 = 1,N 3 = 6) and (N 4 = 2,N 3 = 0). Therefore the 
effects to push the rule table toward chaos is stronger for D4 than D3 of Table HI This point 
will be discussed more quantitatively in the next section. 

At N h =819, 768, 717, 615, 512, 410, 307 and 205 (A = 0.8, 0.75, 0.7, 0.6, 0.5, 0.4, 0.3 
and 0.2), we have carried out replacements experiments for 119, 90, 90, 117, 107, 92, 103 
and 89 rule tables, which are generated randomly using AT/, (A). At all the AT/, points, we 
have succeeded in changing the the patten classes from class III to class II or class I or vice 
versa, by changing the numbers of the group 1 rules. And in many cases, the edge of chaos 
behaviors are observed between them. 
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B. Chaotic and periodic limit of general CA(N,K) 

Let us study the effects of the replacements theoretically in general cellular automata 
CA(N,K). In the general case too, the rule tables are classified into four groups as shown in 
subsection III A. We have denoted these numbers as N(gl), N(g2), N(g3) and N(g4). When 
Nh is fixed, these numbers satisfy the following sum rules, which are the generalization of 
Eq.E3 

N(gl)+N(g2) = K N - 1 , 
N(g3) + N(g4) = K N - 1 (K-l), 
N(g2) + N(g3) = K N - N h , 
N(gl) + N(gA) = N h . 

However the individual number N(gi) suffers from the fluctuations due to randomness. 
They distribute with the mean given by Table IIII1 In this subsection, let us neglect these 
fluctuations. 

The replacements to decrease the number of the group 1 rule while keeping the A fixed 

TABLE III: The classification of CA(N,K) rules into 4 groups. The \x% represent to K-l while h,i 
and I, 1 to K-l. 





rule 


Average number 


group 1 


T(m,H2, -,0, ...,/xjv) = h 


K N ~ l \ 


group 2 


T(m,n2, -,o, ...,hn) = o 


K N -\l-\) 


group 3 


T(m,fj^,...,i,...,fi N ) = 


K N ~\K -1)(1- A) 


group 4 


T(m,n2, —,1, —,fJ>N) = i 


K N -\K-l)\ 



are given by, 

N(gl) ^ N(gl) - 1, N(g2) - N{g2) + 1, 

N(g3)^N(g3)-l, N{gA) ^ N{gA) + 1. (9) 

In CA(5,4), they correspond to RP4 to RP1 of section III A. 

These replacements stop either when N(gl) = or N(g3) = is reached. Therefore 
when N(gl) < N(g3), which corresponds to A < (1 — -L) in A, all the group 1 rules are 
replaced by the group 2 rules. In this limit, quiescent states at time t will never be changed, 
because there is no rule which converts them to other states, while the group 3 rules have 
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a chance to create a new quiescent state in the next time step. Therefore the number of 
quiescent states at time t is a non-decreasing function of t. Then, the pattern class should 
be class I (homogeneous) or class II (periodic), which we call periodic limit. Therefore the 
replacements of Eq. Q push the rule table toward the periodic direction. 

Let us discuss the reverse replacements of Eq. (jSJ). In these replacements, if N(g2) < 
N(g4), which corresponds to < A, all group 2 rules are replaced by the group 1 rules, 
except for the quiescent condition. In this extreme reverse case, almost all the quiescent 
states at time t are converted to other states in next time step, while group 3 rules will 
create them at different places. Then this will most probably develop into chaotic patterns. 
This limit will be called chaotic limit. 

We should like to say that there are possibilities that atypical rule tables and initial 
conditions might generate a periodic patterns even in this limit. But in this article, these 
atypical cases are not discussed. 

Therefore in the region, 

^A<1-1 (10) 

all the rule tables are located somewhere between these two limit, and by the successive 
replacements of Eq. Q and their reverse ones, the changes of the pattern classes take place 
without fail. This is the theoretical foundation of the replacement experiments of previous 
subsection and also explains why in this region the four pattern classes coexist. 

The Eq. © provide us with a method to control the pattern classes at fixed A. The 
details of the replacements of Eq. depend on the models. In the CA(5,4), they have been 
RP4 to RP1 and RC4 to RC1. They will enable us to obtain a rule tables which are difficult 
to generate by the method "random-table method" or "random- walk-through method" and 
lead us to the new understandings on the structure of CA rule tables in section V. 



IV. QUIESCENT STRING DOMINANCE PARAMETER F IN CA(5,4) 

In the previous section, we have found that each rule table is located somewhere between 
chaotic limit and periodic limit, in the region ^ < A < 1 — -4. In order to express the 
position of the rule table quantitatively, we introduce a new quiescent string dominance 
parameter F, which provides us with a new axis (F-axis) orthogonal to A. Minimum of F 
is the periodic limit, while maximum of it corresponds to chaotic limit. In this section, we 
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will determine the parameter F for CA(5,4). 

As a first approximation, the parameter F is taken to be be a function of the numbers of 
the rules D4, D3, D2 and Dl, which have been denoted as A 4 , N 3 , N 2 and Aq, respectively. 
We proceed to determine F(N 4 , N 3 , N 2 , Aq) by applying simplest approximations and as- 
sumptions 

We apply Taylor series expansion for F, and approximate it by the linear terms in N 4 , 
N 3 N 2 and Aq. 

F(N 4} N 3 , N 2 ) ~ c 4 N 4 + c 3 N 3 + c 2 N 2 + Cl N x . (11) 

where c 4 = dF/dN 4 , similar for c 3 , c 2 and c\. They represent the strength of the effects 
of the rules D4, D3, D2 and Dl to push the rule table toward chaotic direction. These 
definitions are symbolic, because N is discrete. 

The measure in the F is still arbitrary. We fix it in the unit where the increase in one 
unit of A 4 results in the change of F in one unit. This corresponds to divide F in Eq. (jlljl 
by C4, and to express it by the ratio c 3 /c 4 (r 3 ), c 2 /c 4 (r 2 ) and c\/c 4 (ri). 

Before we proceed to determine r 3 , r 2 and rq, let us interpret the parameter F geomet- 
rically. Most generally, the rule tables are classified in 1024-dimensional space in CA(5,4). 
The location of the rule table of each pattern classes forms a hyper-domain in this space. 
We map the points in the hyper-domain into 4-dimensional (N 4 , N 3 , N 2 , Ni) space, in which 
they will also be located in some region. We introduce a surface S(N 4 , N 3 , N 2 , Ni) = $ 
in order to line up these points. F-axis is a normal line of this surface. In Eq. (jllj) . we 
approximate it by a hyper plane. 

In order to determine r 3 , r 2 and n, we apply an argument that the class IV rule tables 
are located around the boundary of the class II and class III rule tables. Our strategy to 
determine r 3 , r 2 and rq is to find the regression hyper plane of class IV rule tables on four 
dimensional space, (N 4 ,N 3 , N 2 Aq). It is equivalent to fix the F-axis in such a way that the 
projection of the distribution of class IV rule tables on F-axis, looks as narrow as possible. 
The quality of our approximations and assumptions reflects the width of the distribution of 
class IV rule tables. 

In the least square method, our problem is formulated to find r 3 , r 2 and r 1; which mini- 
mize the quantity, 
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s(r 3 ,r 2 ,r 2 ) = ^^(F^ lassIV (Ni, N^, N%, N[) — F 3 lassIV (N{, N^, N(, N()) 2 ., (12) 



c 4 



where i and j label the class IV rule tables. We solve the equations, dS/dr 3 = 0, dS/dr 2 = 
and dS/dri = 0, which are 

i,3 i,j i,j i,j 

r 3 J2 SNi' j 6N l 2 ' j + r 2 ^(5N^) 2 + nj^ SN^SN? = - $K' J 5N : 



1,3 
2 5 



1,3 h3 1,3 h3 (13) 

r 3 ^ SN^SN^ + r 2 J2 SNpSN? + r x ^(5N^) 2 = - SN^SN^. 

i,3 i,j i,3 i,3 



where SNl' 3 — N\ — Nl, similar for SN^ 3 and 5N% 3 . 

In order to collect class IV rule tables, we have generated rule tables randomly both for 
Nh in the region, 205 < < 819 (0.2 < A < 0.8) and for the numbers of the group 1 rules 
in the ranges, < N A < 6, < N 3 < 33, < N 2 < 72 and < N x < 144. 

This is realized by the two step method. In the first step, we generate rule table randomly 
using the number Nh, which are explained in subsection II B. We should like to notice that 
under this methods, the numbers of the group 1 rules, N^, N 3 , N 2 and N x are distributed 
around, 6A, 33A, 72A and 144A respectively. They are denoted as N*. Then in second 
step, iVjS are determined randomly between zero and their maximum. The iV*'s, which are 
obtained in the first step are changed to their random value by the replacements RC4 to 
RC1 or RP4 to RP1. 

We have generated about 14000 rule tables, and classify them into four pattern classes 
according to their transient length. There are 483 class I, 3169 class II, 10248 class III and 
329 class IV rule tables, respectively. From the 329 Class IV rule tables, the coefficients r« 
are determined by solving the Eq. (|T3*j) . They are summarized in the Table IIV1 where the 
errors are estimated by the Jackknife method. The results show that the coefficients are 
positive, and satisfy the order, 

c 4 > c 3 > c 2 > c\. (14) 

It means that the effects to move the rule table toward chaotic limit on the F-axis are 
stronger for the rules which break longer strings of the quiescent states jl^ . 

The order in Eq. (fPfj) is understood by the following intuitive arguments. If six D4 
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TABLE IV: The optimal and intuitive coefficients r^. 





optimal 


Error 


Intuitive 




0.1563 


0.0013 


0.18182 




0.0506 


0.0007 


0.08333 


n 


0.0195 


0.0002 


0.04167 



rules are included in the rule table, the string of the quiescent states with length 5 will 
not develop. Similarly, if 33 D3 rules are present in the rule table, there is no chance 
that the length 4 string of the quiescent states could be made. These are roughly similar 
situations for the formation of pattern classes. Thus the strength of the D3 rules r 3 
will be roughly equal to 6/33 of that of D4 rules, similar for the strength of the D2 and 
Dl rules. These coefficients are also shown in the Table II VI We call the F parameter 
with these fj's as intuitive F parameter and those determined by solving the Eq. (|13|) as 
optimal one. It is found that the differences between them are not so large. The classifi- 
cation of the rule tables with intuitive F parameter will be discussed in the subsection VII B. 

V. DISTRIBUTION OF THE RULE TABLES IN (A,F) PLANE 

Using the optimal F parameter we plot the position of the rule tables of each pattern 
classes in (A,F) plane. They are shown in Fig. El 

The Figs. 01 shows that the class III rule tables are located in the larger F region; about 
F > 4, while class I, class II rule tables, in the smaller F region; about F < 9, and the 
class IV rule tables are found in the overlap region of class II and class III rule tables; about 
4 < F < 9. These results support the chaotic limit and periodic limit discussed in the 
subsection III B, and shows that at least in the 0.2 < A < 0.75 range, all four pattern classes 
coexist. 

These distributions of rule table in F, are almost independent of A. It means that the CA 
pattern classes are not classified by A, contrary to the results of Ref. jjj, but are classified 
rather well by the quiescent string dominance parameter F. 

Let us discuss the reason why Langton had obtained his results. If the rule tables are 
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FIG. 3: Distribution of rule tables in (A,F) plane. Fig. |3^a): The distribution of 483 class I rule 
tables. Fig.^b): The distribution of 3169 class II rule tables. Fig.|3^c): The distribution of 10248 
class III rule tables. Fig. Efe): The distribution of 329 class IV rule tables. 

generated using the probability A by the "random-table method" or "random- walk-through 



method" 



31 



the numbers of the group 1 rule tables A^4, V3, N2 and N\ are also controlled 
by the probability A. They distribute around of section IV. Then the F parameters are 
also distributed around, 

F x = (6 + 33r 3 + 72r 2 + 144ri)A. (15) 

Therefore in these methods, A and F are strongly correlated. The probabilities to obtain 
the rule tables, which are far apart from the line given by Eq. (|15J1 are very small. When 
A is small, the rule table with small F are mainly generated, which are class I and class II 
CAs. On the other hand in the large A region, the rule tables with large A are dominantly 
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generated, which are class III CAs. The line of Eq. ()15|) crosses the location of class IV 
pattern classes around 0.4 < A < 0.55. This might be a reason Langton has obtained his 
results. But the distribution of rule tables in all (A,F) plane, show the global structure of 
the CA rule table space as in Fig. |3] and lead us to deeper understanding of the structure of 
CA rule tables. 

We should like to stress again that four pattern classes do coexist in a rather wide range 
in A, which are rather well classified by the parameter F not by A. 

VI. CLASSIFICATION OF RULE TABLES IN (A,F) PLANE 

In the Fig. |3J it is found that rule tables are not separated by sharp boundaries. And 
they seem to have some probability distributions. We denote the probability densities of 
class I, class II, class III and class IV pattern classes as P 1 , P 11 , P IH and P IV , respectively 
and proceed to classify the rule tables by using them. 

The equilibrium points F I J^ III {\) of class II and class III rule tables are defined by the 
point where the relation P H (F) = P IH (F) is satisfied. The region in (X,F) plane where P n 
and P I!I coexist in a similar ratio is defined as transition region. The upper points of the 
transition region Fy~ IU , are defined by the points, p n =E^lL and similarly for the lower 
points of the transition region F^ 1-111 , where P 11 and P IH are interchanged. By these three 
points, F^~ H1 , Fl/~ in and F/ 7-777 , we define the phase boundary of the rule tables. 

The distributions of the rule tables in the Fig. |3] show the qualitative probability distri- 
butions. However in order to study the A dependences of F^ 1-111 , F^~ in and p^ 1-111 more 
quantitatively, we generate a rule tables at fixed As. The A points and numbers of the rule 
tables are shown in Table |3 At each fixed A, we divide the region in F into bin sizes 
of SF = 1, and count the number of rule tables of each classes in these bins. From these 
results, we estimate the probability densities P(Fi), where Fi is a middle point of that bin. 

Let us proceed to the classification of the class II and class III rule tables. 

A. Classification of rule tables in 1 < A < 1 — \ 

The probability distributions of P n (F) and P in (F) at A = 0.3 are shown in Fig. QJa), 
and the determination of the Fj/~ IH , F^ 7-777 and p^~ in are demonstrated in Fig. HJb). 
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TABLE V: The data points, numbers of each pattern classes and the transition parameters of 
CA(5,4). The column "Comp" shows the numbers of rule tables which transmit the information 
of initial states. 





Number of the rule tables 


Class I - II 


Class II - III 


A 


N h 


I 


II 


III 


IV 


Comp 


pi-n 

Li 


f i-ii 

Hi 


F i-n 

U 


pii-m 


pii-in 

T? 

Hi 


pii-in 

± TT 

U 


0.125 


128 


202 


1081 


1052 


42 


16 


0.0 


2.0 


3.2 


3.0 


6.6 


9.0 


0.15 


154 


99 


499 


545 


23 


17 


0.7 


2.0 


3.6 


3.5 


6.5 


9.1 


0.2 


205 


141 


728 


1021 


39 


18 


1.2 


2.0 


3.6 


4.7 


6.3 


8.5 


0.25 


256 


98 


479 


949 


19 


8 


0.0 


1.7 


3.8 


4.4 


5.9 


7.5 


0.3 


307 


83 


364 


878 


13 


6 


1.2 


2.0 


3.9 


4.0 


6.0 


7.3 


0.4 


410 


117 


385 


1169 


23 


7 






3.7 


4.5 


5.7 


7.0 


0.5 


512 


53 


341 


884 


23 


7 








4.8 


6.0 


7.2 


0.6 


615 


38 


488 


1316 


61 


10 








4.8 


6.1 


7.3 


0.7 


717 


4 


333 


974 


89 


7 








4.9 


5.8 


6.8 


0.75 


768 


2 


256 


960 


108 


7 








4.5 


5.4 


6.6 


0.8 


819 


2 


81 


1064 


12 


5 












4.3 



For the other A points of Table W\ Fe, Fy and Fl are determined in the similar way. They 
are summarized in the Table 

As already seen in Fig. El the A dependences of the F^~ IH , F^ 1-111 and ar e 
small, in the region 0.25 < A < 0.75. This is also confirmed by the studies at fixed As as 
shown in in the Table IVl 



B. Classification of rule tables in larger and smaller A region 

The P n and P I!I at A = 0.8 are shown in Fig. |S| It should be noticed that there is 
no region of F where P 11 > P 111 . This means that F^ 1-111 and F^~ in disappears. Only 
Fy~ T is determined. 

In order to understand what has changed at A = 0.8, we have studied the A dependences 
of P 11 and P IH in the region A > 0.7. They are shown in Fig. |3 It is found that P HI 
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Probability Distribution min(P" 1 P l ")/max(P" 1 P l ") 




FIG. 4: The probability distributions P 1 , P n and P IU and determination of Fj/~ IH , F^" 111 
and F{/~ HI at A = 0.3. Fig.E^a): P 1 , P n and P IH distributions. Fig. Bp)): An example of the 
determination of Fg~ HI , F^' 111 and F^ 1 ' 111 . 




FIG. 5: The probability distributions P n and P 111 and the determination of , F L and 

Flj 1 - 111 at A = 0.8. Fig-Efa) shows P n and P HI . In Fig.^b), it is found that F| /_m , F^ 111 
could not be obtained. 



gradually increases as A becomes larger but the increase is quite small, while P decreases 
abruptly between A = 0.75 and A = 0.8. As a result P n becomes less than P IH in all 
F regions. This tendency could already be observed in Fig. 0(b), but it is quantitatively 
confirmed by the studies at fixed As. 
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FIG. 6: The A dependences of P and P in the region, A > 0.7 




FIG. 7: The A dependences of P and P in the region, A < 0.3 



In the smaller A region (A < 0.3), the behavior of P n and P IH are shown in Fig. 
In this case, P 11 is gradually increasing as A decreases but the change is small. On the 
contrary, the decrease in P 111 is larger. As a consequence of these changes the transition 
region of the class II and class III pattern classes spread over wider range in F. These 
results are also summarized in Table IS and shown in Fig. |HJ 

In the same way, the classification of the class I and class II rule tables could be carried 
out. The preliminary results are shown in the column F^~ n , 1 and F[~ u of the 
Table El In this case too, it is seen in Fig. Efa), (b) that density of class I rule tables 
decreases as A increases, while that of class II rule tables stays almost constant in A < 0.75. 
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FIG. 8: The phase diagram of CA(5,4). The minimum of F (FMin) and the maximum of F (Fmo,x), 
which are discussed at the end of subsection VI B, are shown by the dotted lines. 



This feature is more quantitatively confirmed by the studies at fixed As. At A = 0.4, there 
disappears the region of F, where P 1 is larger than P n and F[~ n and F^T 11 could not be 
determined, just as in the same way as P n and P IH distributions at A = 0.8. These results 
are also shown in Table IV1 

However we should like to say that the numbers of the class I rule tables, and those of 
class II rule tables in the region F < 3 are not large. Therefore the results may suffer from 
large statistical fluctuations. We think that the classification of class I and class II rule 
table needs more data to get quantitative conclusions, however the qualitative properties 
will not be changed. 

We proceed to the investigation of the classification of rule tables outside of these A 
region. In A < 0.25 region, not all the group 2 rules could be replaced by the group 1 rules. 
Therefore the maximum numbers of group 1 rules N(gV)Max could not become 256, and it 



22 



decreases to zero as A approaches to zero. Then the maximum of F, (-FWr) also decreases 
to zero toward A = 0. 

Conversely in A > 0.75 region, not all the group 1 rules could be replaced by the group 2 
rules. The minimum of N(gl) and therefore the minimum of F, (-F/vfin) could not becomes 0. 
The line Fuin increases until its maximum at A = 1. In Fig- [HI we have schematically shown 
the Fuax and Fmiu lines with dotted lines. We should like to stress that the dotted line 
should have some width due to fluctuations of N4, A3, N 2 and N\ caused by the randomness. 

VII. DISCUSSIONS AND CONCLUSIONS 
A. Transmission of initial state informations 

The computability of the CA is discussed very precisely mainly for elementary CA 
(CA(3,2)) in series of paper from Santa Fe Institute j4j. In this subsection we discuss on 
the simplest problem of transmission of the initial state information to the later configura- 
tions. 

We have found some examples where class II and class IV patterns appear with similar 
probability by changing the initial configurations randomly. An example is shown in the 
Fig. EJ which is the transmission of initial state information to later configurations and is 
similar to the p = 1/2 problem in the CA(3,2). It is interesting to investigate under what 
condition the changes of the pattern classes are taken place. 

We have focus on the difference of patten classes between class II and class IV, because 
in this case differences of the patterns are obvious. These rule tables are found in the wide 
region in A, 0.125 < A < 0.8. The numbers of the rule tables of this property at fixed As are 
also shown in the column " Comp" of Table 

In addition, there are cases where the difference of patters seems to be realized within 
the same pattern classes. In these cases careful studies are necessary to distinguish the 
difference of these patterns. In this article, we have not studied these cases. 
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FIG. 9: An example of the transmission of the initial state information at A = 0.75. Fig. E(a), and 
Fig.EJb) are generated by the same rule table. Only the difference is the initial states configurations 
which are set randomly. 

B. Classification of rule table by the intuitive F parameter 

The methods to determine the coefficients in Eq. (fTTj) are not unique. In the section 
IV, in order to determine them we have used regression hyper plane of class IV rule tables, 
and in order to obtain 329 class IV rule tables, we have generated totally about 14000 rule 
tables. It is a rather tedious task. However optimal set of has been close to the intuitive 
set of r;. 

In this subsection we study the classification of rule tables by the intuitive F parameter. 
The same analyses as sections V and VI are carried out and as the similar figures are 
obtained in this case too, we will show only the distributions of the rule tables in (A,F) 
plane in Fig. El The Figs. El are very similar to the Figs. EJ Therefore the classification of 
the rule table space are almost same as Fig. |Hl except that the Fmqx changes from 17.6 to 



If intuitive F parameter could successfully classify the rule tables for general CA(N,K) 
it would be very convenient, because it reflect the structure of the CA rules and there is no 
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FIG. 10: The distribution of rule tables by using intuitive F parameter. The numbers of the data 
points are same as in the case of FigEl 



need to gather a lot of class IV rule tables, in order to determine r^s. Whether it is correct 
or not must be concluded after the studies of other CA(N,K) lo) ]. 



C. Conclusions and discussions 

We have started to find the mechanism which distinguishes the pattern classes at same 
As, and have found that it is closely related to the structure of the rule tables. The patten 
classes of the CA are mainly controlled by the numbers of the group 1 rules, which has been 
denoted by N(gl). 

In the CA(N,K),in the region < A < 1 — the maximum of N(gl)Max corresponds to 
a chaotic limit, and its minimum N(gl) = 0, to a periodic limit. Therefore in this A region, 
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we could control the patten classes by changing N(gl) without fail. The method for it is 
the replacements of Eq. Q. Using the replacements, we could study the rule tables which 
are difficult to obtain by the "random-table method" or "random- walk-through method" 
of Ref. 0| . This property could be studied quantitatively by introducing a quiescent string 
dominance parameter F. 

In this article, a quantitative studies are carried out for CA(5,4). In this case, the group 
1 rules are further classified into 5 types as shown in Table HI and the classification of rule 
tables is carried out in (A,F) plane as shown in Fig. |H| It is seen that the A dependences of 
the transition region are very gentle, and rule tables are classified better by the F parameter 
rather than by A. It is interesting whether or not the A dependences of the transition region 
depend on the models. 

In the replacement experiments, we have found the edge of chaos (very long transient 
lengths) behavior in many cases. The examples are shown in Fig. |21 Sometimes they are 
observed in some range in N$ or iV 2 . This indicates that in many cases, the transitions are 
second-order like. But the widths in the ranges of N$ or N2 are different from each other, 
and there are cases where the widths are less than one unit in the replacement of RP2 (first- 
order like). It is very interesting to investigate under what condition the transition becomes 
first-order like or second-order like. The mechanism of the difference in the transitions is an 
open problem and it may be studied by taking into account effects of group 3 and 4 rules. 
In these studies another new parameters might be found and a more quantitative phase 
diagram might be obtained. 

These issues together with finding the points where the transition region crosses Fuax 
and F Min lines (dotted lines) in Fig. |Hl and the nature of the transition at these points will 
be addressed in the forthcoming publications. 
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